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Abstract. We provide a pedagogical overview of defect models of structure 
formation. We first introduce the concept of topological defect, and describe 
how to classify them. We then show how defects might be produced in phase 
transitions in the Early Universe and approach non-pathological scaling so- 
lutions. A very heuristic account of structure formation with defects is then 
provided, following which we introduce the tool box required for high pre- 
cision calculations of CMB and LSS power spectra in these theories. The 
decomposition into scalar vector and tensor modes is reviewed, and then 
we introduce the concept of unequal-time correlator. We use isotropy and 
causality to constrain the form of these correlators. We finally show how 
these correlators may be decomposed into eigenmodes, thereby reducing a 
defect problem to a series of "inflation" problems. We conclude with a short 
description of results in these theories and how they fare against observa- 
tions. We finally describe yet another application of topological defects in 
cosmology: baryogenesis. 



1. Introduction 

Phase transitions are ubiquitous in nature. Typically, as the temperature 
of a system drops below the critical temperature, the system makes a tran- 
sition from a state with greater symmetry to one with less symmetry. In 
general, the state with less symmetry is not unique, but can lie anywhere 
in a so-called vacuum manifold. Depending on the topology of this vacuum 
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manifold, defects will form during the phase transition. If the topology of 
the vacuum manifold admits defects, then these defects will inevitably arise 
during the phase transition unless the dynamics is completely adiabatic. 

In the early Universe the temperature was decreasing very rapidly. On 
length scales larger than the Hubble radius, causality prevents the sys- 
tem from maintaining adiabaticity through interactions, and therefore on 
these scales defects will arise in any cosmological phase transition in which 
they are topologically allowed. The Universe has undergone several phase 
transitions. We are quite confident about those which occurred at lower 
temperatures: the confinement transition at a temperature T ~ 10 2 GeV 
and the electroweak symmetry breaking transition at T ~ 10 3 GeV. Uni- 
fied field theories of fundamental interactions predict the existence of other 
transitions at higher temperatures, e.g. a phase transition at T ~ 10 16 GeV 
in Grand Unified Models, the supersymmetry breaking phase transition in 
supersymmetric models, and various compactification transitions in string 
(and M-) theory. 

Since topological defects carry energy density, they will curve space- 
time and can thus act as the seeds for gravitational accretion (see e.g. 
Refs. [1, 2, 3, 4, 5] for comprehensive reviews). Since inside of topological 
defects the symmetry characteristic of the high temperature phase is un- 
broken, topological defects can interact in various interesting ways with the 
surrounding matter and can have an effect on cosmological issues such as 
baryogenesis (see e.g. [6, 7, 8]), magnetic field generation (see e.g. [9, 10]), 
and ultra-high-energy cosmic ray production [11, 12, 13, 14]. 

In these lectures, we first review the classification of topological defects 
and explain why in models with the appropriate topology, defects will in- 
evitably form during the symmetry breaking phase transition. In Section 
3 we discuss some initial applications of topological defects to cosmology. 
We review the domain wall and monopole problems and explain why the 
cosmic string model yields a promising mechanism for structure formation. 

In the following sections we provide a more technical description of 
how high accuracy calculations of structure formation in defect theories 
are performed. We first describe the details of the scalar, vector and tensor 
decomposition (Section 6). This is an invaluable tool in linear perturbation 
theory. Then in Section 7 we introduce the concept of UETC and show their 
general form, assuming isotropy and scaling, but not energy conservation. 
In Section 8 we show how causality limits further the form of the correlators, 
in the large wavelength limit. These results will be important when checking 
upon the numerics. Then in Section 11 we present the UETCs measured 
for cosmic strings, and highlight some of their features. Their remarkable 
novelty is the dominance of the energy density over any other components 
of the stress energy tensor. This property sets strings apart, resulting in 
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a dominance of scalar modes over vector and tensor modes. We present 
some conclusions on string scenarios of structure formation, and also hybrid 
scenarios combining strings and inflation. 

In the final section we illustrate one application of topological defects 
to cosmology which involves microscopic physics rather than gravitational 
accretion: we discuss the basics of defect-mediated baryogenesis. 

2. Defect Classification and Formation 

According to our current particle physics theories, matter at high energies 
and temperatures must be described in terms of fields. Gauge symmetries 
have proved to be extremely useful in describing the standard model of 
particle physics, according to which at high energies the laws of nature are 
invariant under some non-abelian group of internal symmetry transforma- 
tions G = SU(3) C x SU(2)l x U(1)y, where the first factor is the symmetry 
group of the strong interactions, and the second and third factors form the 
gauge group of the Glashow-Weinberg-Salam theory of electroweak inter- 
actions which is spontaneously broken to the U{\) of electromagnetism at 
a scale of T ~ 10 3 GeV. 

Spontaneous symmetry breaking is induced by an order parameter <p 
taking on a nontrivial expectation value < <p > below a certain temperature 
T c . In some particle physics models, ip is a fundamental scalar field in a 
nontrivial representation of the gauge group G. However, <p could also be a 
fermion condensate, as in the BCS theory of superconductivity. 

The transition taking place at T = T c is a phase transition and T c is 
called the critical temperature. From condensed matter physics it is well 
known that in many cases topological defects form during phase transi- 
tions, particularly if the transition rate is fast on a scale compared to the 
system size. When cooling a metal, defects in the crystal configuration will 
be frozen in; during a temperature quench of 4 He, thin vortex tubes of the 
normal phase are trapped in the superfluid; and analogously in a temper- 
ature quench of a superconductor, flux lines are trapped in a surrounding 
sea of the superconducting Meissner phase. 

In cosmology, the rate at which the phase transition proceeds is given by 
the expansion rate of the Universe. Hence, topological defects will inevitably 
be produced in a cosmological phase transition [15], provided the underlying 
particle physics model allows such defects. 

Topological defects can be point-like (monopoles), string- like (cosmic 
strings) [16] or planar (domain walls), depending on the particle physics 
model. Also of importance are textures [17, 18], point defects in space-time. 
Topological defects represent regions in space with trapped energy density. 
These regions of surplus energy can act as seeds for structure formation. 



4 



JOAO MAGUEIJO AND ROBERT H. BRANDENBERGER 



Consider a single component real scalar field with a typical symmetry 
breaking potential 

V{p) = \\{p?-r?f (1) 

Unless A « 1 there will be no inflation. The phase transition will take 
place on a short time scale r < H~ l , and will lead to correlation regions 
of radius £ < t inside of which p is approximately constant, but outside of 
which (f ranges randomly over the vacuum manifold M, the set of values 
of p which minimizes V{(p) - in our example p = ±77. The correlation 
regions are separated by domain walls, regions in space where p leaves the 
vacuum manifold M. and where, therefore, potential energy is localized. 
Via the usual gravitational force, this energy density can act as a seed for 
structure. 

As mentioned above, topological defects are familiar from solid state 
and condensed matter systems. The analogies between defects in particle 
physics and condensed matter physics are quite deep. Defects form for the 
same reason: the vacuum manifold is topologically nontrivial. The argu- 
ments which say that in a theory which admits defects, such defects will 
inevitably form, are applicable both in cosmology and in condensed matter 
physics. Different, however, is the defect dynamics. The motion of defects in 
condensed matter systems is friction-dominated, whereas the defects in cos- 
mology obey relativistic equations, second order in time derivatives, since 
they come from a relativistic field theory. 

Turning to a classification of topological defects, we consider theories 
with an ra-component order parameter p> and with a potential energy func- 
tion (free energy density) of the form (1) with 

<f 2 = t ■ ( 2 ) 

i=i 

There are various types of local and global topological defects (regions 
of trapped energy density) depending on the number n of components of p. 
The more rigorous mathematical definition refers to the homotopy of Ai. 
The words "local" and "global" refer to whether the symmetry which is 
broken is a gauge or global symmetry. In the case of local symmetries, the 
topological defects have a well defined core outside of which ip contains no 
energy density in spite of nonvanishing gradients Vp: the gauge fields 
can absorb the gradients, i.e., D^p = when d^p 7^ 0, where the covariant 
derivative is defined by 

D lx = 8 lx + ieA lx , (3) 

e being the gauge coupling constant. Global topological defects, however, 
have long range density fields and forces. 



COSMIC DEFECTS 



5 



Table 1 contains a list of topological defects with their topological char- 
acteristic. A "v" marks acceptable theories, a "x" theories which are in 
conflict with observations (for rj ~ 10 16 GeV). 



TABLE 1. Classification of cosmologically allowed 
(v) and forbidden (x) defects. 



Defect name 


n 


Local defect 


Global defect 


Domain wall 


1 


X 


X 


Cosmic string 


2 


V 


V 


Monopole 


3 


X 


V 


Texture 


4 




V 



We now describe examples of domain walls, cosmic strings, monopoles 
and textures, focussing on configurations with maximal symmetry. 

Domain walls arise in theories with a single real order parameter and 
free energy density given by (1). The vacuum manifold of this model consists 
of two points 

M = {±7]} (4) 
and hence has nontrivial zeroth homotopy group: 

n (M) + i . (5) 

To construct a domain wall configuration with planar symmetry (with- 
out loss of generality the y — z plane can be taken to be the plane of 
symmetry), assume that 

ip(x) ~ 7] x ^> i] 1 (6) 

tp(x) ~ —7] X <C —T]^ 1 (7) 

By continuity of ip, there must be an intermediate value of x with (p(x) = 0. 
The set of points with ip = constitute the center of the domain wall. 
Physically, the wall is a thin sheet of trapped energy density. The width w 
of the sheet is given by the balance of potential energy and tension energy. 
Assuming that the spatial gradients are spread out over the thickness w we 
obtain 

W V(0) = wXr] 4 ~ - 7] 2 (8) 



and thus 



(9) 
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A theory with a complex order parameter (n = 2) admits cosmic strings. 
In this case the vacuum manifold of the model is 

M = S 1 , (10) 

which has nonvanishing first homotopy group: 

U 1 (M) = Z^1. (11) 

A cosmic string is a line of trapped energy density which arises whenever 
the field <p{x) circles M along a closed path in space {e.g., along a circle). In 
this case, continuity of ip implies that there must be a point with (p = on 
any disk whose boundary is the closed path. The points on different sheets 
connect up to form a line over-density of field energy, a cosmic string. 

To construct a field configuration with a string along the z axis [19], take 
(f(x) to cover M. along a circle with radius r about the point (x, y) = (0, 0): 

(^(r, i?) — rje^ , r 3> f] ^ • (12) 

This configuration has winding number 1, i.e., it covers M exactly once. 
Maintaining cylindrical symmetry, we can extend (12) to arbitrary r 

<p{r,0) = f(r)e i *, (13) 

where /(0) = and f(r) tends to r] for large r. The width w can again be 
found by balancing potential and tension energy. The result is identical to 
the result (9) for domain walls. 

For local cosmic strings, i.e., strings arising due to the spontaneous 
breaking of a gauge symmetry, the energy density decays exponentially for 
r S> . In this case, the energy fi per unit length of a string is finite and 
depends only on the symmetry breaking scale rj 

[i~if (14) 

(independent of the coupling A). The value of /x is the only free parameter 
in a cosmic string model. 

If the order parameter of the model has three components (n = 3), then 
monopoles result as topological defects. The vacuum manifold is 

M = S 2 (15) 

and has topology given by 

n 2 (M)^i. (16) 

Given a sphere S in space, it is possible that <p takes on values in M. 
everywhere on S, and that it covers M. once. By continuity, there must be 
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a point in space in the interior of S with tp = 0. This is the center of a 
point-like defect, the monopole. 

To construct a spherically symmetric monopole with the origin as its 
center, consider a field configuration <p which defines a map from physical 
space to field space such that all spheres S r in space of radius r » rf 1 
about the origin are mapped onto M: 

ip: S r — ► M (17) 

(r,0,yO — (18) 

This configuration defines a winding number one magnetic monopole. 

Domain walls, cosmic strings and monopoles are examples of topological 
defects. A topological defect has a well-defined core, a region in space with 
if £ M. and hence V(ip) > 0. There is an associated winding number which 
is quantized, i.e., it can take on only integer values. Since the winding 
number can only change continuously, it must be conserved, and hence 
topological defects are stable. Furthermore, topological defects exist for 
theories with global and local symmetries. 

Now, let us consider a theory with a four-component order parameter 
(i.e., n = 4), and a potential given by (1). In this case, the vacuum manifold 
is 

M = S 3 (19) 
and the associated topology is given by 



U 3 (M)^1. (20) 

The corresponding defects are called textures [17, 18]. 

Textures, however, are quite different than the previously discussed 
topological defects. The texture construction will render this manifest. To 
construct a radially symmetric texture, we give a field configuration <p(x) 
which maps physical space onto M.. The origin in space (an arbitrary 
point which will be the center of the texture) is mapped onto the north 
pole N of M.. Spheres in space of radius r surrounding are mapped onto 
spheres in A4 surrounding N, with distance from N increasing as r in- 
creases. In particular, some sphere with radius r c (t) is mapped onto the 
equator sphere of M. The distance r c (t) can be defined as the radius of the 
texture. Inside this sphere, f(x) covers half the vacuum manifold. Finally, 
the sphere at infinity is mapped onto the south pole of A4. The configura- 
tion <f(x) can be parameterized by [18] 

¥(x,y,z) = \ cosx(r), sinx(r)-, sinx(r)-, sinx(r)- (21) 
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in terms of a function %(r) with %(0) = and x(oo) = n. Note that at all 
points in space, ip(x) lies in M. There is no defect core. All the energy is 
in spatial gradient (and possibly kinetic) terms. 

In a cosmological context, there is infinite energy available in an infinite 
space. Hence, it is not necessary that x( r ) —> n as r —> 00 • We can have 

X(r) Xmax < vr as r ^ 00 . (22) 

In this case, only a fraction 

Xmax sm (2Xmax) fnn\ 

nw = — (23) 

of the vacuum manifold is covered: the winding number nw is not quantized. 
This is a reflection of the fact that whereas topologically nontrivial maps 
from S 3 to S 3 exist, all maps from R 3 to S 3 can be deformed to the trivial 
map. 

Textures in R 3 are unstable. For the configuration described above, the 
instability means that r c (t) — > as t increases: the texture collapses. When 
r c (t) is microscopical, there will be sufficient energy inside the core to cause 
(p(0) to leave Ai, pass through and equilibrate at x(0) = ir: the texture 
unwinds. 

A further difference compared to topological defects: textures are rele- 
vant only for theories with global symmetry. Since all the energy is in spatial 
gradients, for a local theory the gauge fields can re-orient themselves such 
as to cancel the energy: 

D^ = 0. (24) 

Therefore, it is reasonable to regard textures as an example of a new 
class of defects, semi-topological defects. In contrast to topological defects, 
there is no core, and <p(x) £ M for all x. In particular, there is no potential 
energy. In addition, the winding number is not quantized, and hence the 
defects are unstable. Finally, they exist as long-lived coherent configurations 
only in theories with a global internal symmetry. 

The Kibble mechanism [15] ensures that in theories which admit topo- 
logical or semi-topological defects, such defects will be produced during a 
phase transition in the very early Universe. At high temperatures T 3> T c , 
the symmetry is unbroken. The ground state of the finite temperature ef- 
fective potential [20, 21] (see e.g. [22] for an introductory review) is the 
symmetric configuration 99 = 0. Once the Universe cools below the temper- 
ature T c , the symmetry is broken, and <p(x) rolls into the zero temperature 
vacuum manifold M.. However, by causality there can be no correlation 
between the specific points in M. which are taken on on scales larger than 
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the correlation length £(t). In a relativistic theory the causality bound on 
£is 

m<t, (25) 

where t is the causal horizon. This leads to the formation of defects with a 
mean separation of 

The correlation length £(t) can be determined by equating the free en- 
ergy gained by symmetry breaking (a volume effect) with the gradient en- 
ergy lost (a surface effect). As expected, £(T) diverges at T c . Very close to 
T c , the thermal energy T is larger than the volume energy gain E corr in 
a correlation volume. Hence, these domains are unstable to thermal fluc- 
tuations. As T decreases, the thermal energy decreases more rapidly than 
E corr . Below the Ginsburg temperature Tq, there is insufficient thermal 
energy to excite a correlation volume into the state ip = 0. Domains of size 

£(t G ) ~ A-V 1 (26) 

freeze out [15, 23]. The boundaries between these domains become topolog- 
ical defects. An improved version of this argument has been given by Zurek 
[24]. 

We conclude that in a theory in which a symmetry breaking phase 
transition satisfies the topological criteria for the existence of a given type 
of defect, a network of such defects will form during the phase transition 
and will freeze out at the Ginsburg temperature. The correlation length is 
initially given by (26), if the field (p is in thermal equilibrium before the 
transition. Independent of this last assumption, the causality bound implies 
that £(t G ) < t G . 

For times t > to the evolution of the network of defects may be com- 
plicated (as for cosmic strings) or trivial (as for textures). In any case (see 
the caveats of Refs. [25, 26]), the causality bound persists at late times 
and states that even at late times, the mean separation and length scale of 
defects is bounded by £(t) < t. 

Applied to cosmic strings, the Kibble mechanism implies that at the 
time of the phase transition, a network of cosmic strings with typical step 
length £(tc) will form. According to numerical simulations [27], about 80% 
of the initial energy is in infinite strings (strings with curvature radius larger 
than the Hubble radius) and 20% in closed loops. 

The evolution of the cosmic string network for t > tc is complicated. 
The key processes are loop production by intersections of infinite strings 
and loop shrinking by gravitational radiation. These two processes combine 
to create a mechanism by which the infinite string network loses energy (and 
length as measured in comoving coordinates) . The dynamics of a string in 
an expanding Universe changes as the curvature radius increases. If the 
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curvature radius is smaller than the Hubble radius H~ l (t), the dynamics 
is dominated by the acceleration term and the string will oscillate at rela- 
tivistic speeds. However, if the curvature radius exceeds the Hubble radius, 
the strings are frozen in by the Hubble damping term in the string equa- 
tion of motion. Therefore, if the string correlation length is smaller than 
the Hubble radius, string intersections will be frequent, the network will 
loose a lot of energy to loops and will increase faster than t. However, 
if at any time £(t) exceeds the Hubble radius, the strings will be effectively 
frozen in comoving coordinates and hence in a standard radiation or matter 
dominated cosmology the Hubble radius will catch up to £(t). According to 
this argument [2], at sufficiently late times the correlation leng! th of the 
string network will always be proportional to its causality limit 

Z(t)~t- (27) 

Hence, the energy density Poo(i) hi long strings is a fixed fraction of the 
background energy density p c (t) 

Pooit) ~ l£(t)~ 2 ~ Vt~ 2 (28) 



or 

P^-G^i. (29) 
Pc{t) * y ' 

We conclude that the cosmic string network approaches a "scaling solu- 
tion" in which the statistical properties of the network are time independent 
if all distances are scaled to the Hubble radius. 

Applied to textures (see e.g. [28] for a review), the Kibble mechanism 
implies that on all scales r > to, field configurations with winding num- 
ber nw > n cr are frozen in with a probability p(n cr ) per volume r 3 . The 
critical winding number n cr is defined as the winding number above which 
field configurations collapse and below which they expand. Only collapsing 
configurations form clumps of energy which can accrete matter. For spheri- 
cally symmetric textures, the critical winding n cr is slightly larger than 0.5 
[29]. 

For t > to, any configuration on scale ~ t with winding number nw > 
n cr begins to collapse (if the scale of the texture is larger than t, the Hubble 
damping term dominates over the spatial gradient forces, and the field 
configuration is frozen in comoving coordinates). After unwinding, <p(x) is 
homogeneous inside the horizon. 

The texture model thus also leads to a scaling solution: at all times 
t > to there is the same probability that a texture configuration of scale 
t will enter the horizon, become dynamical and collapse on a typical time 
scale t. 
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3. Defects and Structure Formation: An Overview 

Topological defects are regions in space with trapped energy density. By 
Newtonian gravity, these defects can act as seeds about which the matter 
in the Universe clusters, and hence they play a very important role in 
cosmology. 

As indicated in Table 1, theories with domain walls or with local monopoles 
are ruled out, and those with only local textures do not give rise to a struc- 
ture formation model. As mentioned earlier, theories with domain walls are 
ruled out since a single wall stretching across the present Universe would 
overdose it. Local monopoles are also problematic since they do not interact 
and come to dominate the energy density of the Universe. Local textures 
do not exist as coherent structures with nonvanishing gradient energy since 
the gauge fields can always compensate scalar field gradients. 

Let us demonstrate explicitly why stable domain walls are a cosmologi- 
cal disaster [30] . If domain walls form during a phase transition in the early 
Universe, it follows by causality (see however the caveats of Refs. [25, 26]) 
that even today there will be at least one wall per Hubble volume. As- 
suming one wall per Hubble volume, the energy density pdw of matter in 
domain walls is 

PDw(t) ~ rfr 1 , (30) 
whereas the critical density p c is 

pc = H2 ^G~ m * et ~ 2 - (31) 
Hence, for 7] ~ 10 16 GeV the ratio evaluated at the present time to is 

^(<b)~ ( -5- V (Vto) ~ 10 52 . (32) 
Pc \m p £ J 

The above argument depends in an essential way on the dimension of 
the defect. One cosmic string per Hubble volume leads to an energy density 
p cs in string 

Pes ~ V 2 t~ 2 ■ (33) 

As we have seen above, the cosmic string network approaches a scaling dis- 
tribution (27). Hence, cosmic strings do not lead to cosmological problems. 
On the contrary, since for GUT models with r\ ~ 10 16 GeV 

(JL.) 2 ^ 1Q -6 (34) 
Pc \m p e I 



cosmic strings in such theories could provide the seed perturbations respon- 
sible for structure formation. 
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Theories with local monopoles are ruled out on cosmological grounds 
[31] (see again the caveats of Refs. [25, 26]) for rather different reasons. 
Since there are no long range forces between local monopoles, their number 
density in comoving coordinates does not decrease. Since their contribution 
to the energy density scales as a~ 3 (t), they will come to dominate the mass 
of the Universe, provided rj is sufficiently large. 

Theories with global monopoles [32, 33] are not ruled out, since there 
are long range forces between monopoles which lead to a "scaling solution" 
with a fixed number of monopoles per Hubble volume. 

Let us now briefly discuss some of the basic mechanisms of the cosmic 
string scenario of structure formation. The implementation of the mecha- 
nisms in concrete models will be covered in the following sections. 

The starting point of the structure formation scenario in the cosmic 
string theory is the scaling solution for the cosmic string network, according 
to which at all times t (in particular at t eq , the time when perturbations 
can start to grow) there will be a few long strings crossing each Hubble 
volume, plus a distribution of loops of radius R<^.t. 

The cosmic string model admits three mechanisms for structure forma- 
tion: loops, filaments, and wakes. Cosmic string loops have the same time 
averaged field as a point source with mass [34] 

M(R) = (3Rfi, (35) 

R being the loop radius and (3 ~ 2ir. Hence, loops will be seeds for spherical 
accretion of dust and radiation. 

For loops with R < t eq , growth of perturbations in a model dominated 
by cold dark matter starts at t eq . Hence, the mass at the present time will 
be 

M(R,t ) = z(t eq )f3Rfi. (36) 

In the original cosmic string model [16, 35] it was assumed that loops 
dominate over wakes. In this case, the theory could be normalized (i.e., [i 
could be determined) by demanding that loops with the mean separation 
of clusters d c i accrete the correct mass, i.e., that 

M(R(d d ),t ) = W u M Q . (37) 

This condition yields [35] 

H ~ 10 32 GeV 2 (38) 

Thus, if cosmic strings are to be relevant for structure formation, they must 
arise due to a symmetry breaking at an energy scale i] ~ 10 16 GeV. This 
scale happens to be the scale of unification (GUT) of weak, strong and 
electromagnetic interactions. It is tantalizing to speculate that cosmology 
is telling us that there indeed was new physics at the GUT scale. 
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The second mechanism involves long strings moving with relativistic 
speed in their normal plane which give rise to velocity perturbations in 
their wake [36]. The mechanism is illustrated in Fig. 1: space normal to the 
string is a cone with deficit angle [37] 

a = 8ttGh . (39) 

If the string is moving with normal velocity v through a bath of dark matter, 
a velocity perturbation 

Sv = AirGfivj (40) 

[with 7 = (1 — w 2 )" 1 / 2 ] towards the plane behind the string results. At 
times after t eg , this induces planar over-densities, the most prominent (i.e., 
thickest at the present time) and numerous of which were created at t eq , the 
time of equal matter and radiation [38, 39, 40]. The corresponding planar 
dimensions are (in comoving coordinates) 

t eq z(t eq ) x t eq z{t eq )v ~ (40 x 40v) Mpc 2 . (41) 

The thickness d of these wakes can be calculated using the Zel'dovich 
approximation [40]. The result is 

d ~ G/J,vj(v)z(t eq ) 2 t eq ~ Av Mpc . (42) 

Wakes arise if there is little small scale structure on the string. In this 
case, the string tension equals the mass density, the string moves at rela- 
tivistic speeds, and there is no local gravitational attraction towards the 
string. 

In contrast, if there is small scale structure on strings, then the string 
tension T is smaller [41] than the mass per unit length fi and the metric of 
a string in z direction becomes [42] 

ds 2 = (1 + h 00 )(dt 2 - dz 2 - dr 2 - (1 - 8Gfi)r 2 dy 2 ) (43) 

with 

h 00 = 4G(n — T) In — , (44) 

ro being the string width. Since hoo does not vanish, there is a gravitational 
force towards the string which gives rise to cylindrical accretion, thus pro- 
ducing filaments. 

As is evident from the last term in the metric (43), space perpendicular 
to the string remains conical, with deficit angle given by (39). However, 
since the string is no longer relativistic, the transverse velocities v of the 
string network are expected to be smaller, and hence the induced wakes 
will be shorter and thinner. 
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Which of the mechanisms - filaments or wakes - dominates is deter- 
mined by the competition between the velocity induced by hoo and the 
velocity perturbation of the wake. The total velocity is [42] 

u = _ 27rG ^~ T ) _ 4„Gniry(v) , (45) 
wy(v) 

the first term giving filaments, the second producing wakes. Hence, for small 
v the former will dominate, for large v the latter. 

By the same argument as for wakes, the most numerous and prominent 
filaments will have the distinguished scale 

teqZ{t eq ) X df X df (46) 

where df can be calculated using the Zel'dovich approximation [43]. 

The cosmic string model predicts a scale-invariant spectrum of density 
perturbations, exactly like inflationary Universe models but for a rather 
different reason. Consider the r.m.s. mass fluctuations on a length scale 
2irk~ 1 at the time tu{k) when this scale enters the Hubble radius -ff _1 (i). 
From the cosmic string scaling solution it follows that a fixed (i.e., tn(k) 
independent) number v of strings of length of the order tn{k) contribute 
to the mass excess 5M(k, tu{k)). Thus 

Note that the above argument predicting a scale invariant spectrum will 
hold for all topological defect models which have a scaling solution, in 
particular also for global monopoles and textures. 

The amplitude of the r.m.s. mass fluctuations (equivalently: of the 
power spectrum) can be used to normalize G\x. Since today on galaxy clus- 
ter scales 

^(Mo)~l, (48) 
the growth rate of fluctuations linear in a(t) yields 

^ (k, t eq ) ~ 10- 4 , (49) 
and therefore, using v ~ 10, 

G\i ~ 10" 5 . (50) 



In contrast to the situation in inflationary Universe models, hot dark 
matter (HDM) is not from the outset ruled out as a dark matter candidate. 



COSMIC DEFECTS 



15 



As non-adiabatic seeds, cosmic string loops survive free streaming and can 
generate nonlinear structures on galactic scales, as discussed in detail in 
Refs. [44, 45]. Accretion of hot dark matter by a string wake was studied in 
Ref. [40]. In this case, nonlinear perturbations develop only late. At some 
time t n i, all scales up to a distance q max from the wake center go nonlinear. 
Here 

<7max ~ Gnv-i{v)z(t eq ) 2 t eq ~ Av Mpc , (51) 

is the comoving thickness of the wake at t n [. Demanding that t n \ corre- 
sponds to a redshift greater than 1 leads to the constraint 

G / u>5-10" 7 . (52) 

Note that in a cosmic string and hot dark matter model, wakes form non- 
linear structures only very recently. Accretion onto loops [46] and onto fila- 
ments [47] provide two mechanisms which may lead to high redshift objects 
such as quasars and high redshift galaxies. 

4. Introduction to high precision calculations with defects 

The last decade has witnessed unprecedented progress in mapping the cos- 
mic microwave background (CMB) temperature anisotropy and the large 
scale structure (LSS) of the Universe. The prospect of fast improving data 
has forced theorists to new standards of precision in computing observable 
quantities. The new standards have been met in theories based on cosmic 
innation[48, 49]. Topological defect scenarios [3, 4] have been more challeng- 
ing. However, recently there have been a number of computational break- 
throughs in defect theories, partly related to improvements in computer 
technology. Most strikingly, the method described in [50] showed how one 
could glean from defect simulations all the information required to compute 
accurately CMB and LSS power spectra. In this method the simulations 
are used uniquely for evaluating the two point functions (known as unequal 
time correlators, or UETCs) of the defects' stress-energy tensor. UETCs are 
all that is required for computing CMB and LSS power spectra. Further- 
more, they are constrained by requirements of self-similarity (or scaling) 
and causality, which enable us to radically extend the dynamical range of 
simulations, a fact central to the success of the method. 

This method was applied to theories based on global symmetries. In 
recent work [51, 52] we have shown how the same method could be applied 
to local cosmic strings (see also [53, 54]). In the next few Sections we shall 
describe in detail the simulation and measurement of UETCs which led to 
the work in [51, 52]; as well as present the analytical tools used in this enter- 
prise. The formalism we had to use is unfortunately more complicated than 
[50]. Local strings have an extra complication over global defects, which 
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stems from the fact that we are unable to simulate the underlying field the- 
ory. Instead, we approximate the true dynamics with line-like relativistic 
strings. This is thought to be reasonable for the large scale properties of 
the stress-energy tensor, but we do not have a good understanding of how 
the string network loses energy in order to maintain scaling. 

This leads to two problems. Firstly one is forced to make assumptions 
about which cosmological fluids pick up this deficit. It is often assumed that 
all the strings' energy and momentum is radiated into gravitational waves, 
approximated by a relativistic fluid. This is by no means certain, and it 
may well be that the energy and momentum is transferred to particles [55], 
and hence to the baryon, photon and CDM components. 

Secondly, it is not enough to find correlators for a reduced number of 
stress energy tensor components (two scalar, two vector, two tensor), and 
then find the others by means of energy conservation. If energy conservation 
can be used then one needs to compute 3 scalar, 1 vector, and 1 tensor 
UETCs. If one is not allowed to make use of energy conservation, one has 
to compute 10 scalar, 3 vector, and 1 tensor UETC. 

In the following sections we first give a qualitative description of the 
technical novelties introduced in defect scenarios. We describe defects as 
active incoherent perturbations. We then describe a set of tools with which 
we can perform high accuracy calculations of structure power spectra in 
these scenarios. 

5. Defects as active, incoherent perturbations 

We first focus on the basic assumptions of inflationary and defect theories 
and isolate the most striking contrasting properties. We define the concepts 
of active and passive perturbations, and of coherent and incoherent pertur- 
bations. In terms of these concepts inflationary perturbations are passive 
coherent perturbations. Defect perturbations are active perturbations more 
or less incoherent depending on the defect [56]. 

5.1. ACTIVE AND PASSIVE PERTURBATIONS, AND THEIR DIFFERENT 
PERCEPTIONS OF CAUSALITY AND SCALING 

The way in which inflationary and defect perturbations come about is rad- 
ically different. Inflationary fluctuations were produced at a remote epoch, 
and were driven far outside the Hubble radius by inflation. The evolution of 
these fluctuations is linear (until gravitational collapse becomes non-linear 
at late times), and we call these fluctuations "passive". Also, because all 
scales observed today have been in causal contact since the onset of infla- 
tion, causality does not strongly constrain the fluctuations which result. In 
contrast, defect fluctuations are continuously seeded by defect evolution, 
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which is a non-linear process. We therefore say these are "active" pertur- 
bations. Also, the constraints imposed by causality on defect formation and 
evolution are much greater than those placed on inflationary perturbations. 

5.1.1. Active and passive scaling 

The notion of scale invariance has different implications in these two types 
of theory. For instance, a scale invariant gauge-invariant potential <!> with 
dimensions L 3 / 2 has a power spectrum 

= <|<£ k | 2 > oc k- 3 

in passive theories (the Harrison-Zel'dovich spectrum). This results from 
the fact that the only variable available is k, and so the only spectrum one 
can write down which has the right dimensions and does not have a scale 
is the Harrison-Zel'dovich spectrum. The situation is different for active 
theories, since time is now a variable. The most general counterpart to the 
Harrison-Zel'dovich spectrum is 

P($) = v 3 F^(kri) (53) 

where F$ is, to begin with, an arbitrary function of x = krj. All other 
variables may be written as a product of a power of rj, ensuring the right 
dimensions, and an arbitrary function of x. Inspecting all equations it can 
be checked that it is possible to do this consistently for all variables. All 
equations respect scaling in the active sense. 

5.1.2. Causality constraints on active perturbations 

Moreover, active perturbations are constrained by causality, in the form 
of integral constraints [57, 58]. These consist of energy and momentum 
conservation laws for fluctuations in an expanding Universe. The integral 
constraints can be used to find the low k behaviour of the power spectrum 
of the perturbations, assuming their causal generation and evolution [59]. 
Typically, it is found that the causal creation and evolution of defects re- 
quires that their energy p s and scalar velocity v s be white noise at low k, 
but that the total energy power spectrum of the fluctuations is required to 
go like k A . To reconcile these two facts one is forced to consider the com- 
pensation. This is an under-density in the matter-radiation energy density 
with a white noise low k tail, correlated with the defect network so as to 
cancel the defects' white-noise tail. When one combines the defects energy 
with the compensation density, one finds that the gravitational potentials 
they generate also have to be white noise at large scales [56] . Typically the 
scaling function F^{krj) will start as a constant and decay as a power law for 
x = krj > x c . The value x c is a sort of coherence wavenumber of the defect. 



18 



JOAO MAGUEIJO AND ROBERT H. BRANDENBERGER 



The larger it is the smaller the defect is. For instance x c ~ 12 for cosmic 
strings (thin, tiny objects), whereas x c « 5.5 for textures (round, fat, big 
things) . Sophisticated work on causality [58] has shed light on how small x c 
may be before violating causality. The limiting lower bound x c « 2.7 has 
been suggested. 

Although we will not here have a chance to dwell on technicalities, 
it should be stated that the rather general discussion presented above is 
enough to determine the general form of the potentials for active perturba- 
tions. This has been here encoded in the single parameter x c . We shall see 
that x c will determine the Doppler peak position for active perturbations. 
Doppler peaks are driven by the gravitational potential, so it should not 
be surprising that the defect length scale propagates into its potential, and 
from that into the position of the Doppler peaks. 

5.2. COHERENT AND INCOHERENT PERTURBATIONS 

Active perturbations may also differ from inflation in the way "chance" 
comes into the theory. Randomness occurs in inflation only when the ini- 
tial conditions are set up. Time evolution is linear and deterministic, and 
may be found by evolving all variables from an initial value equal to the 
square root of their initial variances. By squaring the result one obtains the 
variances of the variables at any time. Formally, this results from unequal 
time correlators of the form 



where a denotes the square root of the power spectrum P. In defect models 
however, randomness may intervene in the time evolution as well as the 
initial conditions. Although deterministic in principle, the defect network 
evolves as a result of a complicated non-linear process. If there is strong 
non-linearity, a given mode will be "driven" by interactions with the other 
modes in a way which will force all different-time correlators to zero on a 
time scale characterized by the "coherence time" 8 c (k,r]). Physically this 
means that one has to perform a new "random" draw after each coherence 
time in order to construct a defect history [56]. The counterpart to (54) for 
incoherent perturbations is 



For \rf — 7]\ = | At/ | > O c (k,rj) we have P($(k, rj), At/) = 0. For At/ = 0, we 
recover the power spectrum P(&(k, r/), 0) = P(&(k,r))). 

We shall label as coherent and incoherent (54) and (55) respectively. 
This feature does not affect the position of the Doppler peaks, but it does 
affect the structure of secondary oscillations. An incoherent potential will 



(d>(k, T/)<D(k', rf)) = S(k - k>(*(fc, v)M*(k, t/)) 



(54) 



($(k, r/)$(k', rf)) = <5(k - k')P($(fc, r/), r/ - rf) . 



(55) 
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drive the CMB oscillator incoherently, and therefore it may happen that 
the secondary oscillations get washed out as a result of incoherence. 

6. Tool 1: Scalar, vector, and tensor decomposition 

Having identified the main qualitative novelties in defecy calculations, we 
now proceed to present the set of tools required for performing high accu- 
racy calculations in these scenarios. We start with the decomposition into 
scalar, vector, and tensor components. Let M j,(x) be the defect stress- 
energy tensor. We may Fourier analyze it 

e^(x) = J d 3 fc9^(k)e ik x (56) 

and decompose its Fourier components as: 

e 00 = P d (57) 

Q 0i = ikiv d + ujf (58) 

@ij = p d Sij + (kkj - ^u)n 5 + 

ifki^-J + kjUY) +Il£ (59) 

with tiujf = 0, fcllY = 0, kmfj = 0, and Uf i = 0. The variables {/, v d , p d , IF 9 } 
are the scalars, {uf, Uf} the vectors, and Ufj the tensors. The decomposi- 
tion can be inverted by means of 

v d = -itiQoi (60) 

<4 = (Si-k%)@ 0j (61) 

and 

P s = \®\ (62) 

u s = |(fc*fc»_^)e y (63) 
u v = 

kkjS lm ) - (kk l Sf + kjk l 5™))Q lm (65) 

although simpler recursive formulae may be written. 

This decomposition and choice of harmonics is slightly different from 
the one used in Hu and White [49], and Kodama and Sasaki [60]. The rea- 
son for this is that in [49] and [60] one assumes that all variables represent 



i^Sl-kiVk^Qlj (64) 
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square roots of power spectra. Isotropy can then be assumed. This pro- 
cedure will indeed lead to the right ensemble average P(k) and Ci if the 
source is coherent (or a decomposition into coherent eigenmodes has been 
performed.) However it is not valid for each realization of an incoherent 
source. 

7. Tool 2: Unequal-time correlators 

The unequal time correlators are defined as 

<G^(k, 7 7 )e* /3 (k, rf)) = C^ a/3 (k, r), rf) (66) 

where Q^ u is the stress energy tensor, k is the wave-vector, and rj and rf 
are any two (conformal) times. The UETCs determine all other 2 point 
functions, most notably CMB and LSS power spectra Ci and P(k). 

All correlators between modes at (k, 77) and (k',r/) will be proportional 
to <5(k — k') due to translational invariance. We shall drop this factor in all 
formulae. The correlators can also be functions of k alone, due to isotropy. 
Since conjugation corresponds to k — ► — k, isotropy implies that the corre- 
lators must be real. Because of incoherence the correlators will be generic 
functions of r) and rj '. 

Furthermore, the form of the S+V+T decomposition fixes further the 
form of the correlators. One can always write down the most general form 
of a correlator, and then contract the result with ki or 5ij, wherever ap- 
propriate, to obtain further conditions. For instance for any scalar S and 
vector Vi we could write 

(SVi) = a{k)k (67) 

By contracting with k\ we would then conclude that a(k) = 0. Proceeding 
in this way we can show that cross correlators involving components of 
different type (S, V, or T) must be zero. Furthermore for any two vectors 
Vi, Wi, one has 

(ViW*) = f vw (k, r), r/)(<% - kkj) (68) 

and so instead of 6 functions we have only one. In order to agree with 
conventions used in the literature we shall instead define this single function 
by means of 

f VW (k, V , V ') = J2(ViW*) (69) 

i 

For correlators involving two tensor quantities one can use similar ar- 
guments to prove the general form 

(UjjUlt) = f nTnT (k, V , r/XSijSu - (SikSfl + 5 u 5 jk ) 

-{5ijk k ki + Skihkj) + (5ikkjk + 5ukjk k + 

Sjkhki + Sjikikk) - hkjkkh) (70) 
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and again we have a single function rather than 21. Again in order to comply 
with conventions in the literature we shall define this single function by 
means of 

/ nTnT (M,V) = £<n^> (71) 

ij 

Hence there should be 10 scalar correlator functions 



j-p d p d j-p d v d j-p d p d j-p d Ii s 
jv d v d jv d p d jv d U s 

jpdpd jpdjjS 



(72) 



3 vector correlators: 



drrV 



••• / 



(73) 



and a single tensor correlator function f uTnT . 

In general these functions are functions of (k, r], ?/), and this is indeed the 
case during the matter radiation transition. However well into the matter 
and radiation epochs there is scaling, and these functions may be written 
as: 

f XY (k,r,,r,') = FXY p T X ' ) (74) 

where XY represents any pair of superscripts considered above, and x = krj 
and x' = krf. The above scaling form results from the dimensional analysis 
argument: [6^(x)] = 1/L 2 , [0^(fc)] = L, [<J(k)] = L 3 , [f XY ] = 1/L. We 
have used units where G = c = 1. 

8. Causality constraints 

From analyticity conditions [61] we can use isotropy and symmetry [61, 62] 
to derive general expansions in k around k\ = 0. This is usually very helpful 
to check numerical results. We find: 



(6oo(A:,r?)6oo(fe,?? / )) 


= X 


(e o(M)e 0l (fc,7/)) 


= ihY 


(eoo(M)©y(M')) 


= V5ij + Whkj 


{eoi(k,r])Q 0j (k,7]')) 


— TS-ij ~\~ XJ kikj 


(@ 0i (k,i])@ jk (k,r)')) 


= i[Qki5jk + R(kj5ik + k^ij 




~\~ S ki k j kfc ] 
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(@ij(k,ri)@ k i(k,r)')) = Adijdki + B(5i k 5ji + SilSjk) 

+C{kikj5kl + k k ki5ij) 
+D(kikkSjl + kikidjk + kjkrfik 
+kjkk5il) + Ekikjkkh (75) 

where A, B etc. are functions independent of k. By applying the scalar, 
vector and tensor decomposition we can then find the specific form of the 
correlators in our formalism to obtain relations between the correlators 
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(76) 

We can now derive k — ► constraints on the correlators. These may help 
to complete the functions F XY near the origin, where they may not be 
accessible from simulations due to low sampling of points on the lattice. 
Specifically, we get ratios between the scalar, vector and tensor anisotropic 
stresses and between the momentum and vorticity, 

ju s u s . ju v u v . y nT n T _ 3-2-4 

f vv :fuLO = 1:2 (77) 

to zeroth order and, 
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to first order. 

9. Tool 3: Simulation determination of the correlators 

Having developped the necessary analytical tools we now describe how to 
apply them to simulations. We consider the work in [51]. The simulation 
starts with a network of strings obtained using the algorithm developed 
in [27] which simulates the symmetry breaking of the underlying field by 
selecting random phases for the field at each point in the lattice. Strings 
are identified on points in the lattice where the phase of the field has a non- 
zero winding number. The strings are described by their position X(u, rj) 
where a is a parameter running along the string and rj is conformal time. 
By imposing the following gauge conditions 



where dots and primes denote derivatives with respect to rj and a, respec- 
tively, the equation of motion takes on the simple form of the wave equation 



and can be discretised on the lattice so that all the possible velocities take 
integer values only [63, 64]. This greatly increases the accuracy and speed of 
string network codes based on this algorithm. The network is then evolved 
using the discretised equations of motion and intercommuting relations. 
To simulate the extraction of energy from the system due to the decay 
of the string loops into gravitational radiation and/or particles, loops of 
a minimum size are excised from the simulation at each time step. This 
also ensures that the network scales with respect to the conformal time r\ 
which enables us to extend the dynamical range of the resulting correlation 
functions beyond the limited range covered in the simulation. Fig. 2 shows 
how the correlation length defined by £ = \fjjjpi scales as a function of 
conformal time. Here, p\ is the density of 'long' strings. 

We performed simulations with box sizes ranging from 128 3 to 600 3 , 
with a cut-off on the loop size of two links. Realisation averages were carried 
out with 256 3 boxes once it was determined that the general form of the 
correlators scaled very accurately with box size. To evaluate the UETCs 
from the simulations we selected times in the range O.liV < t < N/4, 
where N is the box size, when we were sure that the string network was 
scaling, and when boundary effects are still excluded by causality. At each 
of the time steps the vectors X and X from each point along the strings in 
the network were stored. From this we obtained the time evolution of the 
network's stress-energy tensor6 Ml/ at each point on the lattice using 



X.X = 



X 2 + X 2 = 1 



(79) 



X 2 + X 2 = 



(80) 




(81) 
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dv 




dv 

Figure 1. Sketch of the mechanism by which a long straight cosmic string S moving 
with velocity v in transverse direction through a plasma induces a velocity perturbation 
Av(dv) towards the wake. Shown on the left is the deficit angle a, in the center is a 
sketch of the string moving in the plasma, and on the right is the sketch of how the 
plasma moves towards the wake with velocity dv in the frame in which the string is at 
rest. 
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Figure 2. Scaling of the correlation length £ for various loop cut off sizes. 
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The Fast Fourier Transforms of all 10 independent components of the 
string stress-energy were then decomposed into irreducible scalar, vector 
and tensor modes using eqns in Section (6). This resulted in all the SVT 
components being computed directly from the simulation without making 
assumptions on energy conservation and on the details of energy dissipation 
from the string network. The drawback of obtaining all the components 
directly in such a manner is that the process becomes computationally 
intensive for even modestly sized simulations, e.g. 256 3 , as, in effect, one 
has to deal with « 10 times the number of variables at each point on the 
lattice. 

By cross-correlating the decomposed stress-energy components from a 
central time with those from all the stored time steps the 14 independent 
UETCs f XY (k,r),r)') were computed. 

In Figs. 3-16 we display the forms of the various functions f XY (k, 25, if) 
where k, rj and rj are in lattice units. This is enough to infer the scaling 
functions F XY (x, x'). One can see that these correlators fall off very quickly 
away from the diagonal, a phenomenon known as incoherence [56] as we 
explained above. Incoherence determines whether or not we have enough 
dynamical range to compute the UETCs: if we see the fall off completly 
clearly we have enough dynamical range! 

A striking feature of our results is the dominance of Ooo over all other 
components (see Figs. 17 and 18). The string anisotropic stresses are in the 
predicted [61] ratios |0 5 | 2 : [O^ 2 : |0 T | 2 of 3 : 2 : 4, as kr -> 0. However 
| Ooo 1 2 ^ I© 5 ! 2 ) and so scalars dominate over vectors and tensors. Also the 
energy density power spectrum rises from a white noise tail at kr ~ into 
a peak at kr « 20, after which it falls off. Sub-horizon modes are therefore 
of great importance. These features consistently appeared for all box sizes, 
and are independent of the cutoff size imposed on the loops. 

10. Tool 4: Decomposition into eigenmodes 

The UETCs c^a^kr, kr') may be diagonalised [50] and written as 

i 

where are eigenvalues. In general, defects are incoherent sources for 
perturbations [56], which means that this matrix does not factorize into 
the product of two vectors v^ v (kT)v a p(kT'). Standard codes solving for 
CMB and LSS power spectra assume coherence. However we see that an 
incoherent source may be represented as an incoherent sum of coherent 
sources. We may therefore feed each eigenmode into standard codes [65] to 
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Figure 3. The function (|p d | 2 ) (k, 25, rf). 




Figure 4. The function (\v d \ 2 )(k, 25, 77). 
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Figure 5. The function (\p d \ 2 )(k, 25, r?). 




Figure 6. The function (|II s | 2 )(fc, 25, rf). 
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Figure 7. The function £\ (\uf \ 2 )(k, 25, rj). 




Figure 8. The function £\ <|nf| 2 )(fc, 25, 77). 
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Figure 9. The function £\ . <|0^ | 2 > (A;, 25, rf) 




Figure 10. The function (p d v d *)(k, 25, 77). 



JOAO MAGUEIJO AND ROBERT H. BRANDENBERGER 




Figure 11. The function {p d p d *){k, 25, 77). 




Figure 12. The function (p d Tl s *)(k, 25, 77). 
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Figure 13. The function (p d v d *)(k, 25, rf). 




Figure 11 The function (p d U s *)(k, 25, r?) 



JOAO MAGUEIJO AND ROBERT H. BRANDENBERGER 




Figure 16. The function £\ (wfuY*)(k, 25, 77). 
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X=kT) 



Figure 1 7. This plot shows the scaling autocorrelation of the energy density, the scalar 
vector and tensor anisotropic stresses. We see that the energy density dominates all the 
other components. 

find the cjfi and PW (k) associated with each mode. The series 

and A^PW(A;) provide convergent approximations to the power spectra. 

In Fig. 19 we show the eigenmodes corresponding to the energy density 
Goo- We see that the leading eigenmode carries the mark of the peak in 
the energy power spectrum at x « 20. The other eigenmodes become very 
small very quickly and change sign with higher and higher frequency. If 
we compute the Ce power spectrum for each of these modes, and sum the 
series, we find quick convergence shortly after 20 modes have been included 
(see Fig. 20). 

11. Some results 

The methods described above supplied a variety of interesting results for 
cosmic strings. They show that local strings, unlike global defects have a 
Doppler peak. The main uncertainty seems to be in the details of how 
strings loose energy (see e.g. [66, 67]). String decay products are clearly 
the most uncertain aspect of cosmic string theory. By measuring the full 
14 UETCs associated with long strings, we assume nothing about decay 
products when extracting information from simulations The simulations 
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Figure 19. Decomposition into eigenmodes of the energy density. 
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will then also place constraints upon the decay products. 

In Fig. 21 we plot y/[£(£ + l)Q/27r], setting the Hubble constant to 
Hq = 50 Km sec -1 Mpc -1 , the baryon fraction to = 0.05, and assuming 
a flat geometry, no cosmological constant, 3 massless neutrinos, standard 
recombination, and cold dark matter. We also superimpose current experi- 
mental points. The most interesting feature is the presence of a reasonably 
high Doppler peak at £ = 400 — 600, following a pronouncedly tilted large 
angle plateau (cf. [68]). This feature sets local strings apart from global 
defects. It puts them in a better shape to face the current data. 

The CMB power spectrum is relatively insensitive to the equation of 
state of the extra fluid. We have plotted results for w x = 1/3,0.1,0.01. 
Dumping some energy into CDM has negligible effect. Small dumps into 
baryon and radiation fluids, on the contrary, boost the Doppler peak very 
strongly. We plotted the effect of dumping 5% of the energy into the radi- 
ation fluid. 

The LSS power spectra on the other hand is strongly dependent on 
w x . In Fig. 22 we plotted the CDM power spectrum P(k) together with 
experimental points as in [69]. The normalization has been fixed by COBE 
data points. We see that the peak of the spectrum is always at smaller 
scales than standard CDM predictions, or observations. However the overall 
normalization of the spectrum increases considerably as w x decreases. 

The CDM rms fluctuation in 8 fa _1 Mpc spheres is ag = 0.4,0.6,1.8 
for w x = 1/3,0.1,0.01. Hence relativistic decay products match well the 
observed a$ « 0.5. On the other hand in 100 /i _1 Mpc spheres one requires 
bias &ioo = fiooV^ioo = 4.9, 3.7, 1.6 to match observations. 

Energy dumps into radiation have no effect on the CDM power spec- 
trum. However if there is energy transfer into CDM or baryons, even with 
w x = 1/3, the CDM power spectrum is highly enhanced. This is due to 
the addition of small scale entropy fluctuations to the usual fluctuations 
gravitationally induced by the stings. We plot the result of a 5% transfer 
into CDM and a 20% transfer into baryons (with w x = 1/3) for which 
6ioo = 2.0,1.5. 

Hence in our calculations local strings have a bias problem at 100 
fo _1 Mpc, although its magnitude is not as great as found in [54]. It de- 
pends sensitively on the decay products, being reduced if the strings have 
a channel into non-relativistic particles, or if there is some energy transfer 
into the baryon and CDM fluid. The main problem with strings in an Q = 1, 
f&6 = 0.05, = CDM Universe is that the shape of P(k) never seems 
to match observations. This may not be the case with other cosmological 
parameters [70, 71]. 
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12. Strings and inflation 

One way of improving upon the previous situation is to consider mixed 
scenarios: strings and inflation. Recent developments in inflation model 
building, based on super symmetry, have produced compelling models in 
which strings are produced at the end of inflation. In such models the cos- 
mological perturbations are seeded both by the defects and by the quantum 
fluctuations. 

A major drawback of inflationary theories is that they are far-removed 
from particle physics models. Attempts to improve on this state of affairs 
have been made recently, resorting to supersymmetry [72, 73, 74, 75, 76, 
77, 78]. In these models one identifies flat directions in the potentials, which 
are enforced by a (super) symmetry. Such flat directions produce "slow-roll 
inflation". In order to stop inflation one must tilt the potential, allowing 
for the fields to roll down. In so-called D-term supersymmetric inflationary 
scenarios, inflation stops with a symmetry-breaking phase transition, at 
which a U(l) symmetry is spontaneously broken, leading to the formation 
of cosmic strings. This is only the most natural of a whole class of models 
of so-called hybrid inflation. Hence a network of cosmic strings is formed 
at the end of inflation. 

In Figs. (23) and (24) we present power spectra in CMB and CDM pro- 
duced by a sCDM scenario, by cosmic strings, and by strings plus inflation. 
We have assumed the traditional choice of parameters, setting the Hubble 
parameter Hq = 50 km sec" 1 Mpc" 1 , the baryon fraction to Oj, = 0.05, 
and assumed a flat geometry, no cosmological constant, 3 massless neutri- 
nos, standard recombination, and cold dark matter. The inflationary per- 
turbations have a Harrison-Zel'dovich or scale invariant spectrum, and the 
amount of gravitational radiation (tensor modes) produced during inflation 
is assumed to be negligible. 

We now summarise the results. 

— The CMB power spectrum shape in these models is highly exotic. The 
inflationary contribution is close to being Harrison-Zel'dovich. Hence 
it produces a flat small £ CMB spectrum. The admixture of strings, 
however, imparts a tilt. Depending on i?si (the ratio of large angle 
anisotropy due to strings and inflation) one may tune the CMB plateau 
tilt between 1 and about 1.4, without invoking primordial tilt and 
inflation produced gravity waves. 

— The proverbial inflationary Doppler peaks are transfigured in these 
scenarios into a thick Doppler bump, covering the region £ = 200 — 600. 
The height of the peak is similar for sCDM and strings, with standard 
cosmological parameters. The Doppler bump is modulated by small 
undulations, which cannot truly be called secondary peaks. By tuning 
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Rsi one may achieve any degree of secondary oscillation softening. 
This provides a major loophole in the argument linking inflation with 
secondary oscillations in the CMB power spectrum [79, 56]. If these 
oscillations were not observed, inflation could still survive, in the form 
of the models discussed above. 
— In these scenarios the LSS of the Universe is almost all produced by 
inflationary fluctuations. However COBE scale CMB anisotropics are 
due to both strings and inflation. Therefore COBE normalized CDM 
fluctuations are reduced by a factor (1 + Rsi) in strings plus inflation 
scenarios. This is equivalent to multiplying the sCDM bias by \/l + i?si 
on all scales, except the smallest, where the string contribution may be 
non negligible. Given that sCDM scenarios produce too much structure 
on small scales (too many clusters) this is a desirable feature. 

Overall "Strings plus inflation" are interesting first of all as an inflation- 
ary model. Its "flat potential" is not the result of a finely tuned coupling 
constant, but the result of a symmetry. Hence in some sense these models 
achieve inflation without fine tuning. The only free parameters are the num- 
ber of inflationary e-foldings, and the scale of symmetry breaking. These 
parameters also fix the absolute (and therefore relative) normalizations of 
string and inflationary fluctuations. 

The combination of these two scenarios smoothes the hard edges of 
either separate component, leaving a much better fit to LSS and CMB 
power spectra. We illustrated this point in this review, but left out a couple 
of issues currently under investigation which we now summarise. 

The CDM power spectrum in these scenarios has a break at very small 
scales, when string produced CDM fluctuations become dominant over in- 
flationary ones. This aspect was particularly emphasized in [82], and there 
is some observational evidence in favour of such a break. An immediate 
implication of this result is that it is easier to form structure at high red- 
shifts [80, 81]. In [83] it is shown that even with Hot Dark Matter, these 
scenarios produce enough damped Lyman-a systems, to account for the 
recent high-redshift observations. 

Another issue currently under investigation is the timing of structure 
formation [80]. Active models drive fluctuations at all times, and therefore 
produce a time-dependence in P(k) different from passive models. The ef- 
fect is subtle, but works so as to slow down structure formation. Hence for 
the same normalization nowadays there is more structure at high redshifts 
in string scenarios [43, 46, 47]. 
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13. Defects and Baryogenesis 

Baryogenesis forms another overlap area between defects in particle physics 
and cosmology The goal is to explain the observed asymmetry between 
matter and antimatter in the Universe. In particular, the objective is to be 
able to deduce the observed value of the net baryon to entropy ratio at the 
present time 

^(t ) ~ 10- 10 (83) 

starting from initial conditions in the very early Universe when this ratio 
vanishes. Here, Aub is the net baryon number density and s the entropy 
density 

As pointed out by Sakharov [84], three basic criteria must be satisfied 
in order to have a chance at explaining the data: 

1. The theory describing the microphysics must contain baryon number 
violating processes. 

2. These processes must be C and CP violating. 

3. The baryon number violating processes must occur out of thermal 
equilibrium. 

As was discovered in the 1970's [85], all three criteria can be satisfied 
in GUT theories. In these models, baryon number violating processes are 
mediated by superheavy Higgs and gauge particles. The baryon number 
violation is visible in the Lagrangian, and occurs in perturbation theory 
(and is therefore in principle easy to calculate). In addition to standard 
model CP violation, there are typically many new sources of CP violation 
in the GUT sector. The third Sakharov condition can also be realized: After 
the GUT symmetry-breaking phase transition, the superheavy particles 
may fall out of thermal equilibrium. The out-of-equilibrium decay of these 
particles can thus generate a nonvanishing baryon to entropy ratio. 

The magnitude of the predicted tib/s depends on the asymmetry e per 
decay, on the coupling constant A of the ns violating processes, and on 
the ratio nx/s of the number density nx of superheavy Higgs and gauge 
particles to the number density of photons, evaluated at the time when 
the baryon number violating processes fall out of thermal equilibrium, and 
assuming that this time occurs after the phase transition. The quantity e is 
proportional to the CP-violation parameter in the model. In a GUT theory, 
this CP violation parameter can be large (order 1), whereas in the standard 
electroweak theory it is given by the CP violating phases in the CKM mass 
matrix and is very small. As shown in [85] it is easily possible to construct 
models which give the right ub/s ratio after the GUT phase transition (for 
recent reviews of baryogenesis see [86] and [87] ) . 
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The ratio ub/s, however, does not only depend on e, but also on 
nx/s(td). If the temperature Td at the time td is greater than the mass 
rax °f the superheavy particles, then it follows from the thermal history in 
standard cosmology that nx ~ s. However, if Td < mx, then the number 
density of X particles is diluted exponentially in the time interval between 
when T = mx and when T = Td- Thus, the predicted baryon to entropy 
ratio is also exponentially suppressed: 



where g* is the number of spin degrees of freedom in thermal equilibrium 
at the time of the phase transition. In this case, the standard GUT baryo- 
genesis mechanism is ineffective. 

However, topological defects may come to the rescue [6]. As was dis- 
cussed at the beginning of these lecture notes, topological defects will in- 
evitably be produced in the symmetry breaking GUT transition provided 
they are topologically allowed in that symmetry breaking scheme. The topo- 
logical defects provide an alternative mechanism of GUT baryogenesis in 
the following way: Inside of topological defects, the GUT symmetry is re- 
stored. In fact, the defects can be viewed as solitonic configurations of X 
particles. The continuous decay of defects at times after td provides an 
alternative way to generate a nonvanishing baryon to entropy ratio. The 
defects constitute out of equilibrium configurations, and hence their decay 
can produce a nonvanishing ub/s in the same way as the decay of free X 
quanta. 

The way to estimate the ng/s ratio is as follows: The defect scaling 
solution gives the energy density in defects at all times. Taking the time 
derivative of this density, and taking into account the expansion of the Uni- 
verse, we obtain the loss of energy attributed to defect decay. By energetics, 
we can estimate the number of decays of individual quanta which the de- 
fect decay corresponds to. We can then use the usual perturbative results 
to compute the resulting net baryon number. 

Provided that mx < Td, then at the time when the baryon number 
violating processes fall out of equilibrium (when we start generating a non- 
vanishing ns) the energy density in free X quanta is much larger than 
the defect density, and hence the defect-driven baryogenesis mechanism is 
subdominant. However, if mx > Td, then as indicated in (84), the en- 
ergy density in free quanta decays exponentially. In contrast, the density 
in defects only decreases as a power of time, and hence can soon dominate 
baryogenesis. 

One of the most important ingredients in the calculation is the time 
dependence of £(t), the separation between defects. Immediately after the 




ee 



m x /T d 



(84) 
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phase transition at time tf (when the defect network is formed), the sep- 
aration is £(tf) ~ A _1 ry _1 . In the time period immediately following, the 
time period of relevance for baryogenesis, £(t) approaches the Hubble radius 
according to the equation [23] 

m ^/X^) 5/4 . (85) 

Using this result to calculate the defect density, we obtain after some alge- 
bra 

defect ~ 0) (OO) 

s rj s 

where ub/s | o is the unsuppressed value oins/s which can be obtained using 
the standard GUT baryogenesis mechanism. We see from (86) that even for 
low values of T^, the magnitude of ng/s which is obtained via the defect 
mechanism is only suppressed by a power of However, the maximum 
strength of the defect channel is smaller than the maximum strength of the 
usual mechanism by a geometrical suppression factor A 2 which expresses 
the fact that even at the time of defect formation, the defect network only 
occupies a small fraction of the volume of space. 

It has been known for some time that there are baryon number violat- 
ing processes even in the standard electroweak theory. These processes are, 
however, non-perturbative. They are connected with the t'Hooft anomaly 
[88], which in turn is due to the fact that the gauge theory vacuum is 
degenerate, and that the different degenerate vacuum states have differ- 
ent quantum numbers (Chern-Simons numbers). In theories with fermions, 
this implies different baryon number. Configurations such as sphalerons 
[89] which interpolate between two such vacuum states thus correspond to 
baryon number violating processes. 

As pointed out in [90], the anomalous baryon number violating processes 
are in thermal equilibrium above the electroweak symmetry breaking scale. 
Therefore, any net baryon to entropy ratio generated at a higher scale will 
be erased, unless this ratio is protected by an additional quantum number 
such as a nonvanishing B — L which is conserved by electroweak processes. 

However, as first suggested in [91] and discussed in detail in many recent 
papers (see [92] for reviews of the literature), it is possible to regenerate a 
nonvanishing ns/s below the electroweak symmetry breaking scale. Since 
there are ub violating processes and both C and CP violation in the stan- 
dard model, Sakharov's conditions are satisfied provided that one can real- 
ize an out-of-equilibrium state after the phase transition. Standard model 
CP violation is extremely weak. Thus, it appears necessary to add some 
sector with extra CP violation to the standard model in order to obtain an 
appreciable ub/s ratio. A simple possibility which has been invoked often 
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is to add a second Higgs doublet to the theory, with CP violating relative 
phases. 

The standard way to obtain out-of-equilibrium baryon number violat- 
ing processes immediately after the electroweak phase transition is [92] to 
assume that the transition is strongly first order and proceeds by the nu- 
cleation of bubbles (note that these are two assumptions) . 

Bubbles are out-of-equilibrium configurations. Outside of the bubble (in 
the false vacuum), the baryon number violating processes are unsuppressed, 
inside they are exponentially suppressed. In the bubble wall, the Higgs fields 
have a nontrivial profile, and hence (in models with additional CP violation 
in the Higgs sector) there is enhanced CP violation in the bubble wall. In 
order to obtain net baryon production, one may either use fermion scat- 
tering off bubble walls [93] (because of the CP violation in the scattering, 
this generates a lepton asymmetry outside the bubble which converts via 
sphalerons to a baryon asymmetry) or sphaleron processes in the bubble 
wall itself [94, 95]. It has been shown that, using optimistic parameters (in 
particular a large CP violating phase A6cp in the Higgs sector) it is pos- 
sible to generate the observed np/s ratio. The resulting baryon to entropy 
ratio is of the order 

^ -aUaT'^fAecp, (87) 

where ay/ refers to the electroweak interaction strength, g* is the number 
of spin degrees of freedom in thermal equilibrium at the time of the phase 
transition, and m t is the top quark mass. The dependence on the top quark 
mass enters because net baryogenesis only appears at the one-loop level. 

However, analytical and numerical studies show that, for the large Higgs 
masses which are indicated by the current experimental bounds, the elec- 
troweak phase transition will unlikely be sufficiently strongly first order to 
proceed by bubble nucleation. In addition, there are some concerns as to 
whether it will proceed by bubble nucleation at all (see e.g. [96]). 

Once again, topological defects come to the rescue. In models which 
admit defects, such defects will inevitably be produced in a phase transition 
independent of its order. Moving topological defects can play the same role 
in baryogenesis as nucleating bubbles. In the defect core, the electroweak 
symmetry is unbroken and hence sphaleron processes are unsuppressed [97], 
provided that the core is sufficiently thick to contain the sphalerons (in 
recent work [98] it has been shown that this is a rather severe constraint 
on workable string- mediated electroweak baryogenesis mechanisms). In the 
defect walls there is enhanced CP violation for the same reason as in bubble 
walls. Hence, at a given point in space, a nonvanishing baryon number will 
be produced when a topological defect passes by. 
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Defect-mediated electroweak baryogenesis has been worked out in detail 
in [8] (see [7] for previous work) in the case of cosmic strings. The scenario is 
as follows: at a particular point x in space, antibaryons are produced when 
the front side of the defect passes by. While x is in the defect core, partial 
equilibration of ng takes place via sphaleron processes. As the back side 
of the defect passes by, the same number of baryons are produced as the 
number of antibaryons when the front side of the defect passes by. Thus, 
at the end a positive number of baryons are left behind. Obviously, the 
advantage of the defect-mediated baryongenesis scenario is that it does not 
depend on the order and on the detailed dynamics of the electroweak phase 
transition. 

As in the case of defect-mediated GUT baryogenesis, the strength of 
defect-mediated electroweak baryogenesis is suppressed by the ratio SF of 
the volume which is passed by defects divided by the total volume, i.e. 

^~SF^| , (88) 
s s 

where (tib/s)\o is the result of (87) obtained in the bubble nucleation mech- 
anism. 

A big caveat for defect-mediated electroweak baryogenesis is that the 
standard electroweak theory does not admit topological defects. However, 
in a theory with additional physics just above the electroweak scale it is 
possible to obtain defects (see e.g. [99] for some specific models). The closer 
the scale ij of the new physics is to the electroweak scale ijewi the larger 
the volume in defects and the more efficient defect-mediated electroweak 
baryogenesis (however, as pointed out in [98], this effect is counteracted 
by the fact that the defect velocity at T = t\ew decreases as r\ decreases). 
Using the result of (85) for the separation of defects, we obtain (for non- 
superconducting strings) 

SF ~ \{^f 2 v D . (89) 
V 

where vd is the mean defect velocity. Typically [98], the resulting value of 
SF is too small for string- mediated electroweak baryogenesis to be efficient. 

Defect-mediated baryogenesis is much more efficient if the defects are 
domain walls (since by purely geometrical arguments the factor SF will 
be much larger). However, as seen in earlier sections of this review, theo- 
ries with topologically stable domain walls are ruled out because the wall 
network would overdose the Universe. What one would like is a theory in 
which a network of walls forms at some time before the electroweak phase 
transition, remains present until after the transition and then decays. As 
has recently been realized [100], this goal may be achieved in theories with 
embedded defects. 
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Embedded defects (see [101] for a recent review) in a theory with vac- 
uum manifold Ai are exact solutions of the field equations which correspond 
to topological defects with respect to a submanifold Ai' of Ai of strictly 
lower dimension. Embedded defects are unstable in the absence of field fluc- 
tuations. However [100] if plasma effects provide effective masses to some of 
the components of the order parameter of the theory such that the space of 
ground states in the plasma reduces to Ai', then the embedded defects can 
be stabilized in the plasma. Embedded strings which can be stabilized in the 
electric plasma before recombination exist in the standard electroweak the- 
ory. It is interesting to investigate if extension of the minimal electroweak 
theory admit embedded walls. Such walls would then yield ideal candidates 
for defect-mediated electroweak baryogenesis. 

It is important to keep in mind that defects occurring at any scale may 
have important consequences for baryogenesis which should be considered 
when exploring the cosmological implications of the models which admit 
such defects. 
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Figure 21. The CMB power spectra predicted by cosmic strings decaying into loop and 
radiation fluids with w x = 1/3,0.1,0.01,0. We have plotted (£(£ + l)Ci/2n) 1/2 in fiK, 
and superposed several experimental points. The higher curve corresponding to w x = 1/3 
shows what happens if 5% of the energy goes into the radiation fluid. 
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Figure 22. The power spectrum in CDM fluctuations for cosmic strings, with 
w — 0.01, 0.1, 1/3. We plotted also the standard CDM scenario prediction and points 
inferred by Peacock and Dodds from galaxy surveys. The top 2 w x = 1/3 curves corre- 
spond to a 5% transfer into CDM, and a 20% transfer into baryons (top). 
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Figure 23. The CMB power spectra predicted by cosmic strings, sCDM, and by inflation 
and strings with Rsi = 0.25, 0.5, 0.75. The large angle spectrum is always slightly tilted. 
The Doppler peak becomes a thick Doppler bump at I = 200 — 600, modulated by mild 
undulations. 
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Figure 24- The power spectra in CDM fluctuations predicted by cosmic strings, sCDM, 
and by inflation and strings with Rsi — 0.25, 0.5, 0.75. We have also superposed the power 
spectrum as inferred from surveys by Peacock and Dodds. 



